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Abstract 

Wc study a two dimensional system of electrons with Rashba coupling in the constant magnetic 
field B and confining potential. We algebraically diagonalize the corresponding Hamiltonian to end 
up with the solutions of the energy spectrum. In terms of two kinds of operator we construct two 
symmetries and discuss the filling of the shells with electrons for strong and weak B. Subsequently, 
we show that our system is sharing some common features with quantum optics where the exact 
operator solutions for the basics Jaynes-Cummings variables are derived from our results. An 
interesting limit is studied and the corresponding quantum dynamics is recovered. 
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1 Introduction 



The spin-orbit coupling, which couples the electron spin and its orbital motion, has been the subject of 
several theoretical and experimental research [1]. This opens the door to developing a new generation 
of electronic spin (spintronics) and presents a new branch of physics of semiconductors. It puts the 
spin of the electron at the center of interest and exploits the spin-dependent electronic properties of 
magnetic materials and semiconductors. The underlying basis for this new electronics is the intimate 
connection between the electron charge and spin. A crucial implication of this relationship is that 
one can have access to spin through the spin property of the electron orbital in the solid. The link 
between the electron charge and spin is expressed by the spin-orbit interaction in semiconductors. 

Novel spin properties arise from the interplay between Rashba spin splitting [2] and further 
confinement of two dimensional (2D) electrons in quantum wires [31 IH [5l [6], rings [3 [8] or dots 
O Uni [HI [12l [131 [in [El- Spin-orbit coupling has also been shown to affect the statistics of energy 
levels and eigenfunctions as well as current distributions [Tni[I7]. The interplay between spin-orbit cou- 
pling and external magnetic fields was analyzed theoretically using random matrix theory [18]. With 
this respect, Schliemann [19] studied the cyclotron motion and magnetic focusing in semiconductor 
quantum wells with spin-orbit coupling. More precisely, the ballistic motion of electrons in III-V 
semiconductor quantum wells with Rashba spin-orbit coupling in a perpendicular magnetic field was 
investigated. Taking into account the full quantum dynamics of the problem, the modifications of clas- 
sical cyclotron orbits due to spin-orbit interaction was explored and the analogy to Jaynes-Cummings 
was established. 

On the other hand, Jaynes-Cummings model is describing the basic interaction of a two-level 
atom and quantized field, which is also the cornerstone for the treatment of the interaction between 
light and matter in quantum optics [20]. It can be used to explain many quantum phenomena, such 
as the collapses and revivals of the atomic population inversions, squeezing of the quantized field 
and the atom-cavity entanglement. Recent experiments showed that Jaynes-Cummings model can 
be implicated in quantum-state engineering and quantum information processing, e.g. generation of 
Fock states (21] and entangled states [22], and the implementations of quantum logic gates [23], etc. 
Originally, Jaynes-Cummings model is physically implemented with a cavity quantum electrodynamics 
system, see for instance |24] . Certainly, there has been also interest to realize Jaynes-Cummings model 
with other physical systems. A typical system is a cold ion trapped in a Paul trap and driven by 
classical laser beams |251 126] where the interaction between two selected internal electronic levels and 
the external vibrational mode of the ion can be induced. 

Motivated by different investigations cited above and in particular [13 [2Z], we develop our proposal 
to deal with different issues. For this, we consider a 2D system in the presence of an external magnetic 
field B and study the quantum dynamics. But, we include the parabolic potential to confine our 
system and Rashba interaction to make contact with quantum optics. Through the Weyl-Heisenberg 
symmetries, we obtain the solutions of the energy spectrum and construct the algebra ■sm(2) as well as 
sm(1, 1). By considering strong and week B, we show that our system reduces to the Landau problem 
for the first case. By using the Heisenberg picture, we derive two copies of the Jyanes-Cummings 
model oscillating with different frequencies. Finally, we recover the results of without confinement 
case [12] in simply way to conclude that our findings are general and deserve different extensions. 
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The present paper is organized as follows. In section 2 we formulate our problem by setting the 
Hamiltonian and choosing the convenient gauge. In section 3, we introduce a series of annihilation 
and creation operators to diagonalize our Hamiltonian, which serves to determine explicitly the ex- 
act eigenvalues and eigenstates. We construct two symmetries and analyze the system behavior by 
distinguishing the strong and week magnetic field cases in section 4. We establish a link with Jaynes- 
Cummings model and make different comments in section 5. Moreover, to show the relevance of our 
results we study a liming case. Finally, we close by concluding our work and giving some perspective. 



2 Solutions of energy spectrum 

We start by formulating our problem to end up with the appropriate Hamiltonian describing the 
system under consideration. Subsequently, we use the algebraic approach to determine explicitly the 
eigenvalues as well as the eigenstates. These will be used to discuss the possibility to fill the shells 
with electrons when the magnetic field is strong and weak. 

2.1 Hamiltonian formalism 



We consider a system of electrons in the presence of a constant magnetic field B = Biz and con- 
fining potential. By taking into account of the Rashba spin-orbit coupling and Zeeman effect, the 
Hamiltonian for a single electron reads as 

-2 

:i) 



H 



vr 



where tt = p-\- is the conjugate momentum and A is the vector potential. A is the Rashba coupling 
parameter, g is the Lande- factor, /Ub is the Bohr magriGton and. c — (Jy,az) are the Pauh matrices. 
To proceed further, we choose the symmetric Landau gauge 



and write the Hamiltonian ([T]) as 
H 
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The algebraic structure of the above Hamiltonian is easily displayed if we adopt the method of 
separation of Cartesian variables. This process suggests to decompose ([3]) into four parts 



H = Hp + Hr + ^uJcLz + ^gfiBBa^ 



(4) 



such that the free part takes the form 
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and the Rashba coupling in magnetic field is given by 



eB 



(6) 



where ujc 



eB 



is the cyclotron frequency, Lz 

= xpy — ypx is the angular momentum and we have set 



the new frequency a.s uj = ^/u^+^uJq. To close this part, we emphasis that (UD is splitting into two 
independent harmonic oscillator Hamiltonian's supplemented by the angular momentum and Rashba 
spin-orbit coupling added to Zeeman term. This convenient form of the Hamiltonian will help us to 
end up with its diagonalization in the simple way. 

2.2 Solution through Weyl-Heisenberg symmetries 

We introduce the standard machinery and techniques to get the solutions of energy spectrum of the 
Hamiltonian (jj]). Instead of directly using the oscillator annihilation operators 



1 f X Uq 



1 



ilr 



(7) 



we work with two new ones, which are linear superposition of Ux and ay, such as 



ad = -y^ (fflx - lay) , 



1 / 

«9 = {o-x + lay) 



(8) 



where Iq = \J~^ being the magnetic length. Note that, and ag are bosonic operators and satisfy 
the relation commutations 

and other relations vanish, from the above operators, one can obtain the useful identities for the 
conjugate momentum 



h 

TTr = 

2ilo 

as well as for the positions 



h(ad- aU +l2(a„ - a 



7r„ 



h 

2h 



h{ad + al^) -l2{ag + al 



^0 / t t 
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(11) 



where we have set = ^1 + , I2 = — and P = These algebraic structures will play 
a crucial role in solving different issues and more precisely in diagonalizing different Hamiltonian's 
entering in the game. 

We start by writing the Rashba Hamiltonian ([6]) in terms of the annihilation and creation operators 
introduced above. Indeed, we have 

HR = H% + Hi (12) 

where these two parts are given by 



ag 
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In the same way, we can diagonalize the free Hamiltonian and angular momentum to finally end up 
with a new form of the Hamiltonian (jj]). This is 



H 



(14) 



To determine the solutions of energy spectrum of the above problem, we solve the eigenvalue 
equation 



H 



tp2 



E 



1p2 



which gives the eigenvalues 



A2 



where the new frequencies are = ^ (w it Wc)- The corresponding eigenstates read as 
and we show that the amplitudes parameterizing these states are given by 



(15) 
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Having obtained the solutions of the energy spectrum, let us briefly discuss how to recover an 
interesting case from what we generated so far. Indeed, in studying cyclotron motion and magnetic 
focusing in semiconductor quantum wells with spin-orbit coupling, Schliemann [19] introduced the 
Hamiltonian type ([1]) without the confining potential. Therefore, to recover the corresponding solutions 
of energy spectrum, we consider the case = or = in the previous equations. This gives the 
eigenvalues 



'fid 



1 



as well as the eigenstates 



\ng,nd,(7) = u^\ng,nd,'l) + Vnl'^g -l,nd- 1,1) 



(20) 
(21) 



where the amplitudes are given by 



J {hwc + giJ-BB) 



2m\'^hwcnd + T {hwc + gi^'sBY 



2mXhiUcnd + | {huJc + g^J■BBy 



(22) 



(23) 
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We mention that Schliemann [19j used the notation a = hX and the quantum number n = Ud- These 
show clearly that our findings are general and make difference with what obtained in |19j . Certainly, 
this will play a crucial role in the forthcoming analysis. 



2.3 Symmetries and shells 

To make comparison with interesting work [27] dealing with confined 2D system in magnetic field, let 
us introduce two symmetries. These concern the dynamical symmetries su{2) and su{l, 1), which can 
be realized in terms of the shell operators. We start with su{2) where the corresponding generators 
can be realized as 

5+ = a^ag, S- = alttd, = ^ (^'^^ 



which verify the commutation relations 

[S+,S^] = 2S,, [S,,S±] = ±S±. (25) 
These give invariant Casimir operator 

C = 1 (S+S^ + 5_5+) + 52 = ^ - i. (26) 

Therefore, to a fixed value = (n^ + ng)/2 of the operator — ^ there corresponds the {2n + 1)- 
dimensional unitary irreducible representation (UIR) of su{2) in which the operator Sz = ^ assumes 
its spectral values in the range —f^<:'y = {rid — f^g)/'^ ^ 



As far as the second symmetry su(l, 1) is concerned, we consider the generators 



satisfying the relations 



[r+, r_,] = -2ro, [To, r±] = ±t±. (28) 

The Casimir operator then is given by 

V = ]^{T+T^+T.T^)-Ti = -i (^M _ 1^ . (29) 

Similarly to the previous case, to a fixed value r/ = (n^ — ng)/2 + 1/2 > 1/2, with Ud — Ug = a'>Q^ oi 
the operator \ + ^ there corresponds a UIR of su{l, 1) in the discrete series, in which the operator 
To = ^ (if + + -^9 assumes its spectral values in the infinite range r/, r/ + 1, 77 + 2, • • • . Alternatively, 
to a fixed value q = —{ud — ng)/2 + 1/2 > 1/2, with Ug — Ud = —ol > of the operator ^ (— "x + l) 
there corresponds a UIR of su{l, 1) in the discrete series, in which the operator Tq = | (— + l) + A''^ 
assumes its spectral values in the infinite range q, q + 1, q + 2, - ■ ■ . 

Below we will see the importance of the both symmetries introduced above when we analyze two 
interesting cases. With these we underline the system behavior with respect to different limit of the 
magnetic field B. We start our analysis by considering the case where B is strong, which is equivalent 
to the limit oJc ^ wq, and gives the frequencies uj'^ ~ Uc and ~ 0. These tell us that the total 
energy can be approximated by 

En, ^ [rid ± ^) ± \91^bB (30) 
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where the quantum number is = 0, 1 • • • . Now by scahng the energy as 

=F ^gf^BB ~ hujc (^rid =^ = ^"d ("^1) 

one can reahze immediately that our system behaves hke a harmonic oscillator in 2D (Landau problem) 
and for a given there is an infinite degeneracy of the Landau levels. We notice that there are type 
of quantum numbers + ^ and — ^, which means that we have two Hilbert spaces. However, both 
spaces are connected via a linear transformation n = + 1 that allows to move from one space to 
another and vice versa. From symmetry point view, this behavior can be regarded as our system has 
ladder states for a discrete series representations of the algebra su{l, 1) labeled by ^(— a ± 1) where 
a = rifi — Ug < iov a < 0. 

For weak magnetic field, which corresponds to the limit oj^ ojq, we can approximate the frequen- 
cies by LJ^ = u)^ ~ uiQ and therefore write the total energy as 

Eng,na - ^0 {nd + ± 1) it ^gfisB. (32) 

To interpret this result let us rearrange it as follows 

En„n, T ^9I^bB ~ hwo (2A ±l)=ex (33) 

which shows clearly that our system becomes now invariant under the algebra su{2) and therefore each 
Landau level has a degeneracy of order (2A it 1). Note that, here also we have two UIR of dimensions 
(2A + 1) and (2A — 1) for the same algebra where the transition between them can be obtained by 
defining p = A — 1. Furthermore, from (j32p we see that we need 2(Ao ± l)(2Ao ± 1) to fill the shells up 
to the value Aq. 



3 Link with Jaynes-Cummings model 

Very recently, it appeared beautiful connections between different areas of physics. Among them, 
we cite the extraordinary bridge between condensed matter physics and high energy physics through 
the link of graphene with quantum electromagnetic [29]. Also, another connection between massless 
Dirac electrons and quantum optics has been established [30]. These links push and motivate to 
look for bridges and establish contacts between different systems. For this purpose, we make contact 
with another area of physics by showing how our system can be linked to quantum optics through a 
mapping between the corresponding Hamiltonian and Jaynes-Communing model. This may help to 
strength a good knowledge of different aspects of quantum optics. 

3.1 Equivalence between models 

To show the relevance of the results obtained so far, we study the presence and absence of the confining 
potential cases in the Hamiltonian system. For first one, we will show that our Hamiltonian Q is 
formally equivalent to two copies of the Jaynes-Cummings models for atomic transition in a radiation 
field, but oscillating with different frequencies. To do our job, we adopt the same method used by 
Ackerhalt and Rzazewski [28] in analyzing the operator perturbation theory in the Heisenberg picture. 
For second one, we derive the corresponding results in simple way from our findings. 
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To proceed further, we need to rearrange our Hamiltonian in order to deal with each part separately 
and establish the associate link. For this, we start by splitting (jl4p into two parts 



where the first one is 



and the second reads as 



H = Hg + Hd 



1 



Hg = hw alttg + ^ + i?^ + -gfisBaz 



(34) 

(35) 
(36) 



Clearly, these two part are completely diff'erent because they are involving different frequencies and 
therefore different oscillations. 

Let us consider the Hamiltonian (j35p and make contact with Jaynes-Cummings model. In doing 
so, we show that (13511 can be written as 



(37) 



Hg = hUJ Mg - ^ (Jz+ C(W) {a^gO + ttgff'^ 

where the two operators Mg and are give by 



M, 



OJr. 



Ng + a^a- + 

2(w - ojc) 

1 , 

- [a^ ± lay) 



1 



and we have set the constants 7 and C('^) ^ 

7" 
C(u;) 



mhw ( ixic 



For later use, it convenient to define an operator as 



-7 + C(^) WqCJ + agG 



(38) 
(39) 

(40) 
(41) 

(42) 



which verifies the commutation relation [Mg, Cg[ = 0. It tells us that there are two constants of motion 
corresponding to the Hamiltonian p7p . This will help in studying different dynamics of the involved 
operators. 

To study the dynamics related to the Hamiltonian ([37]), we introduce the Heisenberg equation of 
motion for the operators at and . These are 



dt ^ 
dt 



Hg,al 



91 



a 



A straightforward calculation leads 



ih— + huj ] a 
dt 

d_ 

'dt 



ih— + hbj — 27 ) o-"*" = C,{ijj)a^gaz 



(43) 
(44) 

(45) 
(46) 
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Using the relations a~^a~^ = a a =0, (T~^crz = — cr"*" and a az = cr , we show that ()46p can be 
written in terms of the constant of motion Cg as 

ih^ + hw~ + 2Cg^ a+ = Ciuj)al. (47) 

From ()35|) and (|17|) . we derive the same second order differential equation for cT+(t) and This 
is given by 

'4t + '""i (•'"'s + + ^^») ( 4 ) = -^'<-' ( 4 ) ■ (^^* 

We are looking for the quantum dynamics, then we need to solve ()45p and (j47p . One way to do is 
to find solutions of the form 

a+{t) = e'^^'/^s+ + e'^^'/''s- (49) 

at(t) = e*'^«^*/'^/+ + e*'^«"*/\ (50) 

where /3^, and are initial time operators. From the shape of the decomposition (|45p and 
(j47p . one can remark that they have an analogy with the solution of ordinary harmonic oscillator. 
Therefore, the solution will be of the form e^^^^^^ and thus after substitution into ()48p gives a second 
order equation for /3(,. This is 

fil-2 {hu~ + Cg) /3g + hLo~ {hu~ + 2Cg) + C^(a;) = 0. (51) 

By requiring the condition [(3g,Cg] = 0, we show that the corresponding solutions under the decom- 
position forms 

^^ = nw~ + (52) 

where are given by 



a^ = hw" + Cg± ^Cl-Ciuj). (53) 
From the decomposition of ag{t) and (7~^{t), one can notice that the operators constants and 



can be obtained by fixing t = in (j49p and (|50p . These give the relations 



g ' "9 

'+ + r 

Injecting the forms ([l9|) - ([50]) into ([l5]) and ()17|l to end up with 



a+(0) = 4 + (54) 
4(0) = C + C- (55) 



+ «5 ^5 = C(w) (s+ + Sg ) (56) 
agS+ + "9'59 = C(w)(/+ + /^). (57) 



These can be solved to obtain the initial operators as 



— "3 -+ lOoj 
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_ TC(^)cT+(0)±«^4(0) 

— — — -+ • i^yj 



an — an 
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Combining all to end up with the final solutions of ()45p and ()47p . These are 



Alu t 



alit) = ^^|(-C(a;V+(0) + a;at(0))e-^/''+(C(a;V+(0)-a+at(0)^^ 

^9 ^9 

They constitute the exact operator solutions for the basic Jaynes-Cummings variables, which have 
been obtained in [25]. 

As far as the second part (j34p is concerned, we apply the same machinery as before to derive 
similar results. Indeed, by introducing the operator 

Md = Nd + a+a- - ^. . 1 (62) 



and the two constants 



7+ = ] {2fuv+ - gfiBB) (63) 



C'(^) = -Aa/^ 1 + - (64) 



2 V uj 



we write the Hamiltonian (j36p as 

Hd = hw+Md - j+a, + C'{uj) (4cj- + ad(j+) . (65) 
Doing the same job to find the required solutions 



a+{t) = -^_^{(c'(c.)4(0)-a+cT+(0))e-^/'' + (-C'(^)ai(0)+a+a+(0))e-^*/n (66) 
a\{t) = { (-C'{uj)a+{Q) + a-at (0)) e-^/'' + (c'(^)a+(0) - a+4(0)) e^."*/^} (67) 



where different quantities are given by 

a± = Cd±^Jcl-Q'^{u) (68) 
Cd = -l^a, + Q'{u^){ala- + aga^). (69) 

These solutions are showing how to obtain the second copy of the Jaynes-Cummings model from our 
findings. In summary, we conclude that our Hamiltonian is equivalent to two copies of the Jaynes- 
Cummings model but oscillating with different frequencies. 

We close this part by noting that, we can easily obtain the dynamic of the complex position from 
the above solutions. Indeed, using (fTT|) to write 

z{t) = k (4(t) + ag(t)) (70) 



and therefore summing up the adjoint time evolution operator of (j6ip and (|67p to end up with the 
dynamics of z{t). The established link shows clearly that our system is sharing some common features 
with quantum optics. Thus, one can use the present system to handle different issues related to 
Jaynes-Cummings model and vice vera. 
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3.2 Limiting case 

The above results show the analogy to Jaynes-Cummings model and therefore allow us to establish a 
relation with the already published work [25]. Now, we study the case where the confining potential is 
absent, which naturally should lead to the Schliemann results jl9] for Jaynes-Cummings model. These 
will be derived in simple way from our findings to show clearly that our work is general and can be 
extended to deal with different issues. 

To recover the dynamics obtained in [TO] for Jaynes-Cummings model, we start our job by fixing 
the frequency as = in the previous results. This requirement leads to the constraints 

uj~ = 0, h = 2, l2 = (71) 

and therefore according to the dynamical equation ()45p or ()6ip we end up with an operator aj time 
independent 

alit) = 4(0) = at. (72) 
However, ([67|) gives a^(t) time dependent 

= { (c'(^c)^+(0) - r,-4(0)) + (-C'(.;,)a+(0) + r+a^ (q)) e^^*} (73) 

where the roots are 

r± = Cd±(C| + C''(^,))^ (74) 

and all involved functions are now in terms of the cyclotron frequency ujc instead of lo. Returning 
back to the definition of different operators to obtain the time evolution of the position operators in 
the Heisenberg picture. This simply is 

x{t) + iy{t) = lo (4 (t) + og) . (75) 

Since the operators ag is time independent, then at t = we have 

x(0) + iy{0) = hag (76) 

and then after replacing, we find 

x{t) + iy{t) = x(0) + iy(0) + ha\{t). (77) 
Prom (jlOp we can express a)^ in terms of the conjugate momentum as 



t ^0 

which leads to the final form of the complex position 



— (tt^ + iiTy) (78) 



x(t) + iy(t) = x(0)+iy(0)+i — r — ( ' ""^ - i2Xha^ 

rj -r^ \UJc m 



-I 



— _ -A y - i2Xha+ (79) 

rj-r. \ujc m J 
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where the operators valued are given by 

= Cd± [Cj + 2\^mnujc) ^ . (80) 

This nothing but the result obtained by Schliemann [19] in dealing with the same system without 
confinement. Thus, it really shows that our findings are important as well as general in sense that we 
can derive other results. 

4 Conclusion 

We have investigated the basic features of confined two-dimensional system with Rashba spin-orbit 
interaction in the presence of an external magnetic field B. This latter allowed us to end up with a 
confining potential along x and y-directions that has been used to deal with different issues. In par- 
ticular, it has been served to split the corresponding Hamiltonian into two parts. This decomposition 
was useful in sense that different spectrum are obtained and lead to the total solutions of the energy 
spectrum. We have shown that those obtained by Schliemann [19] can be derived in the simple way 
from our solutions. 

Using different operators involved in the Hamiltonian, we have realized two dynamical symmetries 
su{2) and sn(l, 1). These together with the strength of magnetic field B allowed us to discuss the filing 
of the shells with electrons. For strong B, we have concluded that our system behaves like a harmonic 
oscillators in 2D with an infinite degeneracy of the Landau levels. However, for week B our system 
becomes invariant under the algebra su{2) and then each Landau level has an finite degeneracy. 

To make contact with quantum optics, we have elaborated a method based on building two Hamil- 
tonian's from the original one. Indeed, by splitting this later into tow parts we have shown that it is 
possible to recover the Jaynes-Cummings model, which is describing a system with two states. This 
has been done by using the Heisenberg dynamics to find the dynamics of the raising Pauli operator 

and creation operators (aj^, a\). After solving different equations, we have ended up with the exact 
operator solutions for the basic Jaynes-Cummings variables those have been obtained in [28]. To show 
the validity of our results, we have derived those obtained by Schliemann [19] as particular cases. 

The present work can be extended to deal with different issues. For instance, we can use the route 
used by Schliemann [TU] to explicitly study the full quantum dynamics. This is based on expanding 
the initial state of the system in terms of its eigenstates. Another alternative is to use the obtained 
results to study different issues related to graphene and spin Hall effect. 
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